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We propose a novel platform for quantum many body simulations of dipolar spin models using 
current circuit QED technology. Our basic building blocks are 3D Transmon qubits where we 
use the naturally occurring dipolar interactions to realize interacting spin systems. This opens 
the way toward the realization of a broad class of tunable spin models in both two- and one¬ 
dimensional geometries. We illustrate the potential offered by these systems in the context of 
dimerized Majumdar-Ghosh-type phases, archetypical examples of quantum magnetism, showing 
how such phases are robust against disorder and decoherence, and could be observed within state- 
of-the-art experiments. 
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Introduction. - In the present work we propose and 
analyze a novel setup for an analog quantum simula¬ 
tor of quantum magnetism using superconducting qubits. 
The scheme builds on the remarkable recent develop¬ 
ments in Circuit QED CHS] in the context of quan¬ 
tum simulation Emu, and especially the 3D Transmon 
qubit m M- The scheme promises a faithful imple¬ 
mentation of many-body spin-1/2 Hamiltonians involv¬ 
ing tens of qubits using state-of-the-art experimental 
techniques. The central idea behind the present work 
is to exploit the naturally occurring dipolar interactions 
between qubits to engineer the desired spin-spin inter¬ 
actions. In combination with the flexibility offered by 
solid-state setups for realizing arbitrary geometry ar¬ 
rangements, this allows us to design general dipolar spin 
models in ladder and 2D geometries. As we will show, 
our scheme competes favorably with present and envis¬ 
aged quantum simulation setups for magnetism with cold 
atoms and trapped ions [T1HT5] . and enables us to ad¬ 
dress some of the key challenges of quantum simulation 
including equilibrium and non-equilibrium (quench) dy¬ 
namics m- Moreover, we note that exploiting dipolar 
interactions to design dipolar spin models is conceptually 
different, and complementary to the remarkable recent 
experiments with superconducting circuits toward real¬ 
izing the superfluid-Mott insulator transition, based on 
wiring up increasingly complex circuits of superconduct¬ 
ing stripline cavities 17]. 

In our analysis we address two of the key aspects of 
the design of our proposed simulator for quantum mag¬ 
netism. First, we present a feasibility study of state-of- 
the-art experimental setups: this includes a discussion of 
the general mechanism to generate dipolar interactions 
between 3D Transmons, combined with ab initio simu¬ 
lations of the coupling strength in our spin model for 
various geometries. Second, we illustrate how state-of- 
the-art setups, composed of up to a dozen qubits and 
characterized by typical disorder and decoherence rates, 



FIG. 1: Circuit QED implementation of spin models, (a) 
Drawing of one half of a waveguide cavity with two microwave 
couplers. Inside the cavity is a piece of sapphire with mul¬ 
tiple Transmon qubits arranged on a triangular lattice. The 
rotation angles of the qubits are chosen such, that only a few 
of them couple to the fundamental mode of the cavity with 
a predetermined coupling strength, (b) The cavity can be 
loaded with an arbitrary subset of qubits e.g. considering the 
ones in the dot-dashed box in (a) realizes a triangular ladder. 
Here, the bond thickness denotes the interaction strength of 
nearest-neighbor (Ji) and next-nearest-neiglibor (J 2 ) interac¬ 
tions. Additional longer-range contributions are represented 
by dashed lines (only the strongest are shown), (c) Sketch 
of the dimerized phase of the extended XY model in Eq. (JT|). 
The phase can be understood as a solid of local triplet states, 
denoted by the shaded areas. 


are already able to demonstrate paradigmatic signatures 
of quantum magnetism. In particular, we show how a 
dimerized phase |18l , a valence-bond-solid reminiscent of 
the Majumdar-Ghosh state widely discussed in the con¬ 
text of quantum spin chains [19j . can be realized and 
probed with current technologies. 

The system dynamics we are interested in, is described 
by a generalized XY Hamiltonian of the form 

H/h = ^ ■m 1 M {s + s - + hc) + £ hjS , (1) 
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where 5“ are spin operators at the lattice site j, r,j is 
the distance vector between i and j , the last term de¬ 
scribes a disordered transverse field, and the inter-qubit 
couplings Ja = J If' ’| 3 J are in frequency units. The key 
element of our implementation is the realization of dif¬ 
ferent patterns of quantum frustration [181 by tuning the 
form of the interaction couplings (Fig. [2]). As discussed 
below, the latter display a rich dependence as a func¬ 
tion of the dipole angles 6j (see Fig. §>): this depen¬ 
dence is essential to tune the coupling between different 
spins from positive to negative, or to (approximately) 
set it to 0. Even more crucially, the 1 /1 r tJ - 1 3 depen¬ 
dence (see Fig. [2^) can be exploited to further modify 
the magnitudes of the couplings. In the theory of ID 
and ladder systems [201ECO , where the dipolar interaction 
is effectively local, this allows independent tunability of 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
exchange, a fundamental ingredient to realize bond-order 
solids [191 E2 . We focus specifically on this case below, 
showing how, within our proposal, such states of matter 
are robust against both disorder and qubit decoherence. 
The setup can be straightforwardly extended to 2D ge¬ 
ometries, where qualitatively new features emerge due to 
the non-locality of the dipolar couplings: in particular, 
the dynamics of Eq. (IT]) can by-pass the Mermin-Wagner- 
Hohenberg theorem ]23j El] , and thus support phases of 
matter with true long-range order even at finite temper¬ 
ature [25]. We emphasize that all of these phenomena 
are directly accessible within our proposal thanks to the 
naturally occurring dipolar interactions, while such in¬ 
teractions would be challenging to be implemented via 
wiring. 

Model system. - We now describe how the many- 
body dynamics of Eq. (jT|) can be realized in our setup, 
using state-of-the-art circuit QED technology, with the 
design flexibility and long coherence times 1/n (k < 
2-7T x 100 kHz) of Transmon qubits [T21 15] ■ This is 
schematically illustrated in Fig. [I] Several Transmon 
qubits (in red) fabricated on a piece of sapphire (light 
blue area) are mounted inside a waveguide cavity (grey 
box). The Transmon qubits can be fabricated in an es¬ 
sentially arbitrary lattice configuration with locally con¬ 
trollable orientation. The waveguide cavity around the 
qubits is used to readout the state of selected qubits and 
apply a drive to a subset of qubits, providing means for 
both adiabatic state preparation and probing (see next 
section). This selectivity can be again achieved by par¬ 
tially rotating the qubits inside the cavity by a few de¬ 
grees, such that their dipole moment has a finite overlap 
with the electric field in the cavity. To add more flexibil¬ 
ity and mitigate disorder, each qubit can be fabricated 
on an individual piece of sapphire, and can be probed 
at room temperature [26] , For simplicity these individ¬ 
ual pieces are not shown in Fig. [T] From fabrication one 
can expect a disorder of about 3-4% in the Josephson 



FIG. 2: Microscopic finite-element simulation (see text) of the 
dependence of the coupling Ji } j between two qubits. Panel 
(a): distance dependence at 9i = 6j = 90° for three different 
antenna length. The dashed line indicates = 0. The inset 
displays a typical result from HFSS simulations mi, with the 
dipoles surrounded by the computed electric field (intensity 
decreasing from red to blue). Panel (b): angular dependence 
illustrated for the case 9i = 9j = 8 at distance r = 1.5 mm. 
The solid line plots Eq.[2]with Jo,r m extracted from best fits 
in (a). 


energy which one should be able to reduced to less than 
1% (around 27 t x 30 MHz spreading in hj) by preselect¬ 
ing the qubits according to their measured normal state 
resistance. 

Designing the interaction. - The unique feature of 
our approach is the design flexibility in the interqubit in¬ 
teraction which results from the dipole-antenna structure 
of the Transmon. As one can expect, two dipole anten¬ 
nas in the near field interact like two magnetic spins 128] . 
By designing the shape and size of the antenna we can 
realize large interaction strengths (<7^- ss 27 t x 100 MHz) 
at inter-qubit distances of about 1 mm. The interaction 
between the qubits ultimately comes from an effective ca¬ 
pacitance between the antennas, similar to 0®, with an 
angle and distance dependence akin to magnetic dipoles. 

To confirm and quantitatively assess the inter-qubit 
dipolar interactions, we have performed a finite element 
study to determine the coupling strength and dependence 
on distance and angle. These simulations are carried out 
by ab initio numerically solving the 3D Maxwell equa¬ 
tions on finite grids (HFSS software mi), and as such 
provide an extremely accurate quantitative benchmark 
for our modeled inter-qubit couplings. The quantitative 
insights gathered from the simulations are of basic impor¬ 
tance to underpin the effects of the specific configuration 
that would be accessible in realistic setups, where the 
influence of the surrounding waveguide cavity and finite 
size effects of the Transmon antenna have to be system¬ 
atically understood. A set of results can be seen in Fig. [2] 
(see also [25] for additional results). Assuming a dipole- 
dipole interaction with an additional component due to 
the dispersive coupling of both qubits to the cavity with 
strength g (determined from independent simulations), 
we expect the following spatial dependence: 


2 cos 0 — 3 cos 9 
12 “ Jo<i ’" yvyj 


*1 COS 02 (0 s 
r ) 3 ^ ' 

1 m ) 


















3 


The first term takes into account the qubit-qubit inter¬ 
action mediated via the fundamental cavity mode with 
(pi , <p 2 the orientation of the qubits relative to the ori¬ 
entation of the electric field in the cavity. Higher order 
cavity modes contribute with a similar term with a cou¬ 
pling g' depending on the mode structure; as these are 
typically very far detuned, we neglect them in the follow¬ 
ing. The second term is the direct dipole-dipole interac¬ 
tion between the qubits with r = |r 12 | being the distance 
between their centers. Both, the qubit-cavity interaction 
strength (which depends on the location of the qubit in 
the cavity) and the dipole-dipole interaction Jo get mod¬ 
ified by the length of the antenna, giving rise to the term 
d m in the Eq. [2j with d m the normalized antenna length. 
The remaining variables in this equation are the qubit 
cavity detuning A « 2tt x 1.5 GHz and the angles #i,0 2 . 
These are the angles formed by the two dipoles with re¬ 
spect to a line connecting their centers, as can be seen in 
Fig.gb). r m corrects for finite size effects of the dipole 
antennas. 

One can see in Fig. [2ja) that the qubit-qubit interac¬ 
tion strength in our numerical simulations very closely 
follows the analytical expression Eq. [2j The only fitting 
parameters in this expression are Jo, which is the same 
for all three datasets, and the r m , which is adjusted for 
each length. The dependence on the distance for differ¬ 
ent antenna sizes can be seen in Fig. |2](a). The interac¬ 
tion strength between the qubits for a parallel orienta¬ 
tion is negative for short distances and falls of as 1/r 3 . 
For larger distances it becomes positive due to the ad¬ 
ditional coupling mediated off-resonantly via the cavity. 
This leads to the effect that qubits at a given distance 
(in the case of our cavity about 3.5 mm) do not interact 
with each other as these terms exactly cancel out. [44] 

By rotating one qubit around the other, as shown in 
Fig. §b), we can see that the interaction strength goes 
from a negative value for a parallel orientation to a twice 
as large positive value for collinear qubits, as expected 
for a dipole-dipole interaction. The analytical curves 
agree well with the numerical simulations and demon¬ 
strate that the spatial dependence of the interaction be¬ 
haves like a magnetic dipole-dipole interaction. Further¬ 
more we note that the cavity mediated term can be fully 
suppressed by rotating the qubits perpendicular to the 
electric field of the cavity mode, as shown schematically 
in Fig. [I] The finite-element simulations are in very good 
agreement with a simple analytical circuit model (see 
[29]), which already captures the important interaction 
features. 

Frustrated spin ladders and dimer phases. - The 
setup discussed above paves the way to the realization 
of a broad class of lattice spin models JT8]. The key chal¬ 
lenge on the theory side is to identify many-body phe¬ 
nomena which can be realized within this toolbox and, 
crucially, are stable against the most common imperfec¬ 
tion sources such as disorder and decoherence. Here, we 


show how a archetypical state of matter, a valence bond 
solid akin to the Majumdar-Ghosh (MG) state in frus¬ 
trated spin chains Haul EM, can be realized and 
probed with currently available technology. Specifically, 
we are interested in the model in Eq. |l]) on a triangular 
ladder geometry (see Fig. 00- To simplify our discussion, 
since dipolar interactions are effectively local in ID, we 
consider only NN and NNN couplings, that is: 

B\;u\ / h = ^(J 1 ^r +1 + J 2 5+5r +2 + h.c.). (3) 

i 

This model J22[ iili is known to exhibit a phase transition 
from a superfluid phase (Ji J 2 ) to a dimer phase (DP) 
around J 2 /J 1 — 0.33, which is stable up to relatively 
large couplings [3(2. This DP can be easily understood 
at the MG point J 2 /J 1 = 0.5, where the ground state 
factorizes onto a product of local triplets [22] : 

L/2 

\*)0S = G(|t2j-l4-2j)+ |'l'2j —lT 2j )) (4) 

i =1 

as illustrated in Fig. 0b). Beyond being a textbook ex¬ 
ample of a valence bond solid in ID, the DP has triggered 
further interest in connection to the physics of the Peierls 
instability in Hubbard models .32]. In general, the DP is 
indicated by a bond-order parameter (BOP): 

l- 1 

D a = { J2 Df), Df = (-1 ys?s? +1 , a = x,z (5) 
i =1 

which indicates the spontaneous formation of triplets 
along the odd-even bonds on an open chain. Deep into 
the DP, the phase can be identified in an even simpler 
way by monitoring the bond correlation of a single pair 
of spins in the middle of the chain, B z = D z L j 2 (for 
L = 4n,n G N)[45]. In Fig. [3](a), we compare the be¬ 
havior of the dimer correlation B z (red line) against the 
BOP (dashed line) as a function of the coupling ratio 
J 2 /J 1 for a chain of L = 192 sites, obtained using density- 
matrix-renormalization group simulations [33] 146] . Both 
quantities are maximized in the DP, as expected; in ad¬ 
dition, close to the MG point, the difference between the 
BOP and B z is only at the few percent level even for 
very small system sizes of L = 8 spins, providing an ideal 
observable to detect the DP even in relatively small sam¬ 
ples. Measurement techniques to access these observables 
exploiting current technology are discussed in Ref. 129] . 

A crucial point is, to which extent the DP and the 
corresponding signatures on small size systems are stable 
against disorder. Qualitatively, we expect that, close to 
the MG point, the phase is stable against weak disorder 
for hj < J 2 , as breaking a single triplet will cost energy 
~ J 2 . We have performed an exact diagonalization study 
of the resilience against disorder, which is included as 
a Gaussian distribution of hj around 0 with spread Sh 
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FIG. 3: Static properties of Eq. (0). Panel (a): Comparison 
between the BOP D z for a large system (L=192 spins, dashed 
black line), and the bond correlation in the middle of the chain 
for samples of different sizes. Panels (b): Effects of disorder on 
the dimer phase. Up to disorder Gaussian spreadings of order 
5h ~ 0.45Ji, the bond correlation is weakly unaltered i.e. the 
phase is stable against local perturbations. The blue curve 
represents a conservative estimate for typical experimental 
setups. Here, L = 10 sites, results show averaging of 500 
disorder realizations per point. 


(Ji setting the energy unit). In Fig. 0 d, we show the 
behavior of the bond correlation as a function of J 2 /J 1 for 
a system of L = 10 sites and different disorder strength 
and average values. Up to a spread of around 0.45, the 
signatures of the dimer phase are very stable. Given that 
typical disorder strengths would be Sh ~ 0.2/0.3, this 
makes the observation of the DP realistic with current 
fabrication technology. 

Adiabatic state preparation and effect of dissipation. 
- We finally discuss how the DP can be realized in 
realistic settings via adiabatic state preparation. The 
system is initially in the vacuum state, where no ex¬ 
citations are present, that is, j’P) in = ]/[. ® |X)j0?]. 
An external drive is used to inject excitations, and can 
be modeled by a time-dependent Hamiltonian Hd = 
QjSj + Jfj Aj(t)Sj. The dynamics is well described 
by a master equation in Lindblad form (34] : 

m = -^w T (t),p(t)]+ ( 6 ) 

+ \ E [2 C jP (t)Cj - p(t)cjc 3 - cjc,p(t)\ 

Z 1=1,2 L 

where Hx(t) = H + Hd(t), the jump operators C 7 = 
V K jSf and Cj+L = .Sf describe single-spin de¬ 
cay and dephasing, respectively, with constants «,7 < 
100 kHz. In our ramps, we employed a staggered driving 
field, ilj = ((1 + (—1) J ) + (>yL/ 2 )U, which is provided 
by the local spin addressability using angle tuning (see 
Fig. 0a)), and performs significantly better than other 
choices. The specific shape of the sweep parameters we 
choose is illustrated in the inset of Fig. 0i: first, the de¬ 
tunings A j are slowly switched off while ramping up Q. 
Subsequently, the driving is switched off in order to reach 
the desired final state [45] , 

We quantitatively address the full system dynamics in¬ 
cluding disorder and decoherence by numerically solv¬ 




FIG. 4: Adiabatic state preparation. Panel (a): state prepa¬ 
ration of the dimer phase for a system with L = 6 , 2 J 2 = Ji = 
2n x 100 MHz, and different decay rates (in 27rx kHz units), 
with 7 = k. The upper (lower) panel denotes the bond x 
( z ) correlation averaged over 2000 disorder realizations, with 
Sh/Ji = 0.25. The straight lines denote the ideal case. Up to 
decoherence rates of order 100 kHz, the final observables at 
the end of the ramp are within 2% of the expect one. The in¬ 
set shows how the ramping parameters (in units of Ji ) evolve 
with time. Panel (b): same as in (a) for fixed k = 2n x 100 
kHz, and different values of J 2 (in units of Ji). Outside of 
the DP (red dot-dashed line), the order parameters have very 
different final values, as the final state is not dimerized. 


ing Eq. [ 6 ] for systems up to i = 10 spins. In Fig. 0 

(a) , we plot the average magnetization per site (blue 
dashed) and the bond correlation (both in the x - and 
z-basis, black and red line, respectively) as a function 
of time, with a system including decoherence rates with 
k = 7 = 27t x 10,100, 500 kHz, and a disorder strength 
5h = 2n x 25 MHz. The effect of decoherence along 
the ramp is minimal: the final order parameters devi¬ 
ate only a few percent from the expected value even for 
k = 2tt x 100 kHz, making the phase detectable for typi¬ 
cal experimental conditions. For values of the ratio J 2 /J 1 
outside of the DP, dimerization does not occur (Fig. 0 

(b) ). Quantum optimal control techniques [35] ISS] could 
eventually be employed to further improve the ground 
state preparation. 

Conclusions and outlook. - We have shown how to 
use the naturally occurring dipolar interactions in 3D su¬ 
perconducting circuits to realize a platform for analogue 
quantum simulation of XY spin models. The possibil¬ 
ity of realizing arbitrary lattice geometries with locally- 
tunable dipole moments, in combination with their large 
interaction strength, opens the door to the investigation 
of a series of phenomena in quantum magnetism in both 
ID [37, 38] and 2D [25], complementing the remarkable 
developments in cold atom and trapped ion systems Hu¬ 
ng. Our ideas are not limited to Transmon qubits, but 
could be implemented with, e.g., Xmon qubits fB] or 
Fluxonium qubits coupled to an antenna [39]. It would 
be interesting to explore these developments in view of 
realizing Hamiltonian dynamics for surface code architec¬ 
tures [401 or as a building block for coupled cavity array 
experiments 03 m. 
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In the first part and second part of this appendix the details of numerical finite element simu¬ 
lations and the lumped element circuit model are explained. In the third part the Hamiltonian 
for two capacitively coupled Transmon qubits inside a cavity is derived. More details on the 
qubit state measurement and flux tuning are given in part four. In part five, we present an 
extended discussion of the dimer order parameters utilized in the main text to characterize 
the dimer phase. 


HFSS SIMULATIONS 


In this section, finite element electromagnetic simula¬ 
tions that were performed using the ANSYS HFSS [T] 
software are explained in more detail. Qubits are mod¬ 
eled as linear LC resonators (see also .2]). The capac¬ 
itance of the resonator originates from the capacitance 
between two rectangular films that are separated by dis¬ 
tance d. These pads couple to electromagnetic waves like 
a dipole antenna. The Josephson junction of the qubit 
is modeled by a lumped element inductance connecting 
the two antenna pads. Fig. [5] shows the electric field dis¬ 
tribution on a sapphire substrate that holds the qubit 
structure together with structure variable names used in 
this note. 

It should be noted that finite-element calculations of 
a structure with objects that are orders of magnitude 
different in dimensions requires plenty of computational 
resources. In this work, extremely fine meshing is needed 
in areas where accurate understanding of the physics of 
small objects (like qubits) is needed. Fig. [6] shows how 
different mesh sizes are used for different structure size 
to accurately calculate the eigenmodes of the electromag¬ 



FIG. 5: Electric field vectors of the qubit mode plotted on the 
sapphire substrate with the qubit at the center. A rainbow 
color code is used to show the field intensity (with red being 
high intensity). 


netic fields of two qubits in a cavity. 

In this study, the coupling strength is defined as the 
minimum difference in the mode frequencies of two qubits 
when the inductance of one qubit is swept in order to 
tune its mode frequency in and out of resonances with 
the other(see Fig. 14). 


Single qubit inside cavity 

In order to better understand the effect of the an¬ 
tenna geometry on the qubit properties, a few cases have 
been studied. First, the width of the antenna pads w 
was changed and the corresponding capacitance was ex¬ 
tracted as shown in Fig. [7] It can be concluded that the 
capacitance increases linearly as a function of the pad 
width. 

Then, the coupling strength of the qubit to the cavity 
g was calculated as a function of position of the qubit 
with respect to the side wall of the cavity. The coupling 
strength has a cosine form as expected from a standing 
wave in the cavity. The results are shown in Fig. [8] 



FIG. 6: Mesh used to calculate the electromagnetic fields 
around the qubits. The left half of the image shows the mesh¬ 
ing needed for the cavity that contains a sapphire substrate 
with two qubits on it. The right half shows meshing around 
one of the qubits. The coordinate system used in this paper 
is defined in the lower left. 
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FIG. 7: Extracted capacitance of the qubit as a function of 
the width of the antenna paddles. The dashed line is a linear 
fit to the data. 
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FIG. 8: Coupling strength g of the qubit as a function of its 
position inside the cavity. The qubit was placed in the center 
of the cavity in the y and z direction and was moved along 
the x axis. The dashed line is a cosine fit to the data. 


The coupling strength g was also studied as a function 
of the pad height h as shown in Fig.[9]where the qubit was 
located at the center of the cavity. It can be seen from 
the figure that the coupling strength increases linearly as 
a function of the antenna height. This is in agreement 
with dipole antenna physics. 


Two qubits in a cavity 

Fig. [lO] shows how the coupling strength varies when 
one of the qubits is kept at a constant distance (center- 
to-center) of 1.5 mm from the other while the angle of the 
vector that connects the centres of two qubits is changed. 


FIG. 9: Coupling g of the qubit to the cavity as a function of 
its antenna length. Since a change in antenna length will also 
change the resonance frequency, the width of the antenna has 
been adjusted to keep the resonance frequency constant. The 
dashed line is a linear fit to the data. 
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FIG. 10: Coupling as a function of center-center vector 
angle. The two antennas are parallel and their center to center 
distance was kept at 1.5 mm. The dashed line shows a fit to 
the circuit model (see section Q|. 


In our simulations, the coupling strength is defined to be 
a positive number. However, in this particular case, the 
sign of the coupling strength has been changed to nega¬ 
tive for values below « 40 degrees. This is to highlight 
the fact that the sign of the coupling is different in the 
two regions. This sign change can be seen from HFSS 
simulations which show parallel/anti-parallel field vec¬ 
tors for the lowest mode, indicating a positive/negative 
coupling. As explained in section [J the results are in 
strong agreement with the proposed circuit model. 

In another study, the coupling strength was evaluated 
when moving one qubit with respect to the other along a 
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FIG. 11: Coupling J as a function of distance change along 
x axis. The center to center separation along the Y axis is 
kept at 1 mm and one of the qubits is moved along the x 
axis with respect to the other. The results suggest double 
suppression of coupling in agreement with the circuit model 
(dashed curve). 


FIG. 12: Circuit model corresponding to two qubits in a 
cavity. Values of the four capacitors drawn in red decay as 
1 /r 2 where r is the center to center distance between pads of 
the antenna. This models the dipole-dipole coupling of two 
qubits. The two purple structures on Q1 and Q2 represent 
the physical shapes of the qubits. 


straight line. In Fig. ED the coupling strength is shown 
for a qubit moving along the x axis while the distance 
along y is kept at 1 mm. This geometry is inspired by 
the circuit model prediction that there should be two 
zeros along this line. A fit to the circuit model is also 
shown in the graph. 


CIRCUIT MODEL 


In this section, details of a circuit model that repre¬ 
sents two-qubits inside a cavity are explained. The model 
is used to evaluate the coupling strength between the el¬ 
ements. 


Circuit diagram 


Dissipation-less LC resonators are used to model the 
cavity and qubits (see also j2j). The nonlinearity of the 
qubits is not taken into account as it does not affect the 
coupling strength. Capacitors are used in order to ac¬ 
count for the coupling between elements. The coupling 
between each qubit and the cavity is represented by a 
coupling capacitor C cav . As qubits are assumed to cou¬ 
ple to each other through their dipole antenna, a total 
of four capacitors are used to represent coupling of each 
antenna pad to the two pads of the other qubit (see Fig. 
121 - 


Coupling strength 


In order to calculate the coupling between elements, 
the following method is used: The impedance of the sys¬ 
tem ( Z in ) seen through one of the qubits (here Q i) is 
calculated. The poles of the impedance are the reso¬ 
nance modes of the total system, in this example three. 
By sweeping the inductance L\ (of Q\) we can see the 


avoided crossing between the modes (see Fig. 14). Cou¬ 


pling is defined as the minimum distance between the two 
modes at the avoided crossing. 

A further step was taken by relating the inter-qubit 


coupling capacitances (as in Fig. 12) to the physical dis¬ 
tances among the qubit antenna poles. This distance 
dependence is assumed to be of the following form [3|: 


Cab = (2e/7r) In 


= (2e/7r) In 


(«l — b L )(a R — b R ) 
(or — b L )(dL — b R ) 
(w + s) 2 
(w + s) 2 — w 2 


(7) 


where e is the average over dielectric constants of the 
substrate and the air (or vacuum), w = w a = Wb is the 
width of antenna paddles and s is the separation between 
closest edges as shown in the inset of Fig. [l3l The pa¬ 
rameters dL, d R , b R , b R are defined in Fig. [13] 

As this study concerns square shaped paddles (see top 
right inset of Fig. 13), it is not possible to apply equation 


[7]to any arbitrary relative position of the two qubits with¬ 
out complicated corrections. Therefore, just to demon¬ 
strate the general behavior, a simplified relation (Eq. [8] 
follows Eq.[7]when the separation s becomes much larger 
than the paddle width w - this effect is shown in Fig. P 
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FIG. 13: Capacitance between two co-planar paddles (solid 
line) as a function of the separation s according to equation 
[7] The corresponding structure is shown at the bottom. The 
dashed line shows the quadratic decay used in this study 
to estimate the coupling capacitance. The structure at top 
right shows which capacitors in the circuit model follow the 
quadratic decay. 



FIG. 14: Example for three modes of a system with two qubits 
and a cavity.The minimum distance of two modes is defined as 
the coupling 2 Jij if it is two qubits and as 2g if it is the qubit 
cavity coupling. Note that despite the fact that Q 2 has 9 nH 
of inductance, the crossing happens at lower inductances due 
to the coupling to Q\. Dashed line shows how resonance fre¬ 
quency of a typical LC circuit would behave without coupling 
to the other resonators. 


for the coupling is chosen which assumes that the capac¬ 
itance decays quadratically as a function of the distance: 


Cab — n 


( 8 ) 


where a is a coefficient used to adjust the coupling 
strength. 

This means that the coupling strength Jij can be cal¬ 
culated based on the physical position of the two qubits 
with respect to each other. In this study, Q 1 is assumed 
to be fixed in position and Q 2 changes its distance and 
angle with respect to Q\. 

Assuming that the two qubits are placed parallel with 
respect to each other, we have calculated the coupling 
Jij as a function of the center-to-center distance. In 
Fig. [Ts] the results are shown for two qubits without a 
mediating cavity. As can be seen, the coupling is heavily 
suppressed along a straight line. This line can be thought 
of as when the dipole electric field from Qi is horizontal, 
that is perpendicular to the dipole antenna of Q 2 and 
hence, there will be no coupling. This is typical behavior 
of dipoles and is well expected. However, by adding a 


cavity to the model (see Fig. 16), one can observe that 


the straight line shape of suppression will no more hold 
and instead, a curved line emerges. 

Adding the cavity has certain consequences on how two 
qubits couple. For example, if one moves Q 2 along the 
x-axis while it is held at y = 1 mm (or y = — 1 mm), 
the coupling might be suppressed twice. This is verified 
by HFSS simulations as seen in Fig. ED Our lumped- 
element model shows behavior that is consistent with 
HFSS simulations as already demonstrated in Figs. [lO] 
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FIG. 15: Contour plot showing the coupling strength between 
two qubits Jij as a function of the distance between them, 
in the x and y direction, in absence of the cavity. The color 
scale is logarithmic to highlight the rapid change in Jij. In 
such a case, the dipole physics dominates and the coupling 
will be suppressed along a straight line. The white area is not 
accessible due to the dimensions of the qubits. 
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HAMILTONIAN DERIVATION 


Interacting Transmons 


Let us consider two interacting transmons as depicted 


in Fig. 17 The Lagrangian describing this system is given 
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FIG. 16: Contour plot showing the coupling strength between 
two qubits Jij as a function of the distance between them, 
in the x and y direction, in presence of the cavity. The color 
scale is logarithmic to highlight the rapid change in J;,,. The 
cavity distorts the straight line of suppressed coupling opening 
new possibilities for engineering of the coupling between two 
qubits. The white area is not accessible due to the dimensions 
of the qubits. 



FIG. 17: Lumped element circuit model for two capacitively 
coupled superconducting qubits. 


with d> T = (<f>i, ( f> 2 )- For small phase fluctuations, we 
can expand the cosine and get H = Hq + Hi , with 


where 


H 0 

Hi 


1 


QtC 1 Q + — 1 T 


-£ 


E J $4 


r _ (Ci+ C Q -Cq \ 

\ -c Q C 2 + C Q )' 

r -l = I ( C 2 + Cq Cq \ 
D\ Cq Ci + Cq J’ 

r-i _ (Von 0 \ 

^ 0 \/a 2 ) ’ 


( 11 ) 

( 12 ) 


D ee CiC 2 + CiCq + C 2 C Q and 1/c* = Ef /$§. We 
now quantize this Hamiltonian using the usual ladder 
operator replacement for Qj and which is 

Qj = *<9j PF ( a ] “ a i)> = $ j PF ( a I + %)’ ( 13 ) 


where Qj PF = yJhojjCj S /2 , $| PF = sJhujjOij/2 , = 

, C| ff = . Here, uq, w 2 are the two solutions 

of |£ _1 — w 2 C| = 0. This quantization gives 


H S3 ^ uqrq - ^ Ajn 2 + Ai 2 (a|a 2 + h.c.), (14) 

2 = 1,2 2 = 1,2 


r-i(i) //FxZPF \ 4 

where Aj- = (-^-) , A 12 = ^§~Qf PF Q 2 PF , and 

= a j a 3 . Notice that this Hamiltonian has been writ¬ 
ten in a basis of local modes. Diagonalization of the 
quadratic part via a Bogoliubov transformation gives rise 
to a Hamiltonian of the form [2J 


H « ^2 iirii - ^2 °-i n i - Xi 2 nin 2 . (15) 
2 = 1,2 2 = 1,2 


by 


1= Ek + £l 1 -i , 2 ) 2 + (9) 


+ £^co S (|l 


•hr 


with 4>o = h/(2e). The Hamiltonian, H = JT =1 2 Qi®i~ 
L can then be written as 

H = 1 -QtC~ 1 Q - £ Ef cos f, (10) 

i=l,2 ^ 0/ 


Interacting transmons inside a cavity 


Let us now consider the circuit depicted in Fig. 18 


which is a simplified version of Fig. 12 with the capaci¬ 
tive coupling network between the qubits reduced to one 
capacitor. This lumped element model will provide the 
right Hamiltonian but does not capture the angle and 
distance dependence explicitly. We first write the La- 
grangian of the system, which reads 


+ 2 ^^ ~ 2L + _ 

+ £ \ C o& - < m 2 - 


(16) 
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FIG. 18: Lumped element circuit for two interacting super¬ 
conducting qubits in a cavity resonator. 


The associated Hamiltonian H = Hq + Hi can be written 
as equations © and (fl2|), with 


(Ci + Cg + Co -C Q -Co \ 
C = C Q C 2 + C q + Co -Co 

V - Co -Co C + CoJ 

/1/ii 0 0 \ 

C~ x = 0 1/L 2 0 

\ 0 0 l/L ) 


with the renormalized constants are related to the origi¬ 
nal energies as 


Aj 

= —Aj(l — £7 r), 

(21) 

A» 

= Aj(l + 4e^), 

(22) 

Jij 

\ i ^iR^jR 

- AiJ + 2A i ’ 

(23) 


where i = 1, 2. The Hamiltonian (20) can be readily gen¬ 
eralized for N qubits by letting i = 1,N, and con¬ 
veniently diagonalized into the form written in Eq. (15). 
Also, it can be written in a spin language, by considering 
sufficiently large qubit anharmonicities. This allows us 
to replace the bosonic operators by spin 1/2 operators, 
so the effective Hamiltonian for two capacitively coupled 
superconducting Transmon qubits inside a cavity is given 

by 


H cS J U (StS~ + h.c.). (24) 

i i^Lj 


Note that the angle and distance dependence of the qubit 
coupling is hidden in the Jij in this derivation. 


Here, Li = Quantizing this Hamiltonian gives 


QUBIT STATE MEASUREMENT AND FLUX 
TUNING 


H « ^2 UiUi - ^2 A iU 2 i + (17) 

i=l,2,R i—1,2 

+ E + h.C.), 

with similar expressions for the coupling constants as in 
Eq. @ Here, we will consider qubit and cavity far 
detuned. However, the qubit-cavity coupling may be suf¬ 
ficiently large to mediate virtual transitions between both 
qubits through the cavity. This effect will be relevant 
when ~ A12, where A j = —Uj + lor. Performing 

a canonical transformation 

H —► e s He~ s , (18) 

with 

-S = E £ iFi( a i a R - a R a i), (19) 

i=l,2 

and taking £m ss — ff , we obtain, to second order in the 
transformation, an effective Hamiltonian (in a rotating 
frame with respect to the resonator) 

H e s ~ ^2 ^ iTli ~ E + ( 20 ) 

*=1,2 £=1,2 

+ 'y ( Jij(ajaj + h.c.), 
i&= 1,2 


An advantage of a circuit QED setup is the ability to 
measure the state of a predetermined but otherwise ar¬ 
bitrary subset of qubits in the lattice. A measurement 
of (Sf) can be realised in circuit QED by probing the 
cavity transmission in the limit of a dispersive qubit cav¬ 
ity coupling. With the newly developed quantum limited 
amplifier technology it is nowadays possible to even de¬ 
tect quantum jumps of an individual qubit HE]. Using 
this measurement procedure one has to take into account 
the strong interaction between the qubits, which means 
that a direct measurement of (Sf) in the uncoupled basis 
is not possible [Q. We would rather project onto one of 
the eigenstates of the coupled system which might be use¬ 
ful for an initial evaluation of the interaction but would 
not tell us anything about the actual properties of the 
state. 

This effect can be avoided by effectively switching of 
the coupling by detuning a subset of qubits from each 
other and all other qubits by more than J in a time much 
faster than 1 /J in order to be non adiabatic. This de¬ 
tuning can be achieved by realizing the qubits one wants 
to measure with a squid instead of a single Josepshon 
junction. The flux through these qubits can be adjusted 
with the help of small coils just outside the cavity. Simu¬ 
lations and preliminary experiments show that a current 
that is sufficiently large to generate the desired flux can 
be switched in less then 10 ns. Changing the current in 
one of these coils will effect more than one qubit though 
with different magnitude. Adjusting the currents in the 
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FIG. 19: From top to bottom: Scaling of D z , B z , B zz as a 
function of J 2 /J 1 for different system sizes. 


different coils accordingly using e.g. an optimal control 
algorithm should in principle allow us to generate any 
desired flux at the qubit location even though there are 
interdependencies. 


FIG. 20: Scaling of the B zz order parameter as a function of 
the inverse length of the system. 


the center of the chain B z . As already stated in the text, 
these results further confirm that the bond correlation 
well captures the DP indicating that a dimer is indeed 
formed in the center bond of the chain. 

The local order parameter can be further simplified 
by taking the difference between the center bond an a 
neighbor, that is: 


B zz = D z l/2 - D 


'L/2+V 


(25) 


This quantity has been extensively used as a true order 
parameter in numerical simulations for the bond order 
phase of extended Hubbard models (see, e.g., 0), which 
have directly related to the DP here. In the lower panel 
of Fig. [l9j we illustrate the scaling of B zz as a function 
of L and different couplings. Despite its simplicity, this 
combination of correlation functions, which can be di¬ 
rectly accessed in the setup discussed in the main text, 
can serve as an even better indicator of the DP with re¬ 
spect to B z . This is further demonstrated by the results 


in Fig. 20 there, we show how the finite-size scaling of 
B zz is also able to approximately detect the critical point 
between the DP and the superfluid phase in our model. 
An additional discussion on these types of order param¬ 
eters can be found in Ref. [5] 


SCALING OF THE BOND ORDER PARAMETER 
AND BOND CORRELATIONS 

We provide here some additional information on the 
bond order parameter and correlation used in the main 
text to characterize the dimer phase (DP). The main 
message here is that, while the DP is rigorously char¬ 
acterized by the true order parameters D a , simpler cor¬ 
relation functions, which can immediately be accessed in 
the experiments, already qualitatively demonstrate the 
existence of the DP. 

As a first figure of merit, we show in top two panels of 
Fig. [19] a comparison for different system sizes of the true 
bond order parameter D z , and the bond correlation in 
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